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Трансформации силь тяжести. Задача восстановления потенциала силь тяжести в боль- 
ших областях возникла еще в трудах Клеро о фигурах равновесия в виде сфероидов, которьте 
обосновали гравиметрию как науку. Последующие работь Стокса, Пуанкаре и 
Молоденского указали путь ее решения в виде смешанной граничной задачи для уравнения 
Лапласа. Йз-за приближенности граничньх условий зта задача гарантирует 
удовлетворительную точность вьтчислений потенциала только в локальньїх областях «малой 
мерьм» |Якимчик, 20011. 

Другой путь вьтчислений потенциала основан на том, что, согласно |Юнг, 1936), аномалий 
силь тяжести считают в первом приближениий значениями «вертикальной» производной зтой 
аномалий, т. е. реальнье значения силь тяжести заменень гармонической функцией У... Зто 


позволило на базисе теории гармонических функций возвести всю теоретическую надстрой- 
ку современной гравиметрий. Но погрешность такой заменьі, как указано в |Якимчик, 2001| 
и ниже, растет с размерами области сосредоточения аномальньгх масе. Вследствие сказанно- 
го для надежного восстановления потенциала при региональньх исследованиях следуєт во- 
спользоваться иньми методами, например ЇДубовенко, 2009). 

Гравиметрическая сьемка определяет набор значений модуля градиента потенциала и си- 


льт тяжести (МГПСТ) 5Ї,  ргад и, не м(5 ) на некоторой замкнутой дневной поверхности 5, 
которая целиком охватьваєт тяготеющие массьт в области С поиска трансформант. Прямая 
его трансформация не имеет отношения к линейньтм граничньм задачам математической фи- 
зики и является нетривиальной задачей из-за неизвестного направления вектора 8 з: При 
приближенном решений зтой задачи ее граничнье даннье можно задавать на незамкнутой, 
но довольно большой поверхности | Алексидзе, 1965. 


Трансформацию значений гравитационного поля определяют путем решения граничной 
задачи непосредственно для нелинейного дифференциального оператора А (2) сильт тяжести 


в области задания 5 з ІЧерньй, 1970): 


А(а) п. да(х, у- а"(я,х; уз (ж, ) з м(х; ), хєеС,і-і3, (0) 
3 3 

где а" (з, й ) -- У со8 а А со8 а з 293 ра - геометрический фактор поля тяжести, 
ічі уч іч ХОХ; 8 


хх, ) -- -Алу(вгай сх, ) -- заданная функция распределения масс, а, -- угол между направ- 
лением вектора 8), и і-й координатной оси. Зто уравнение учитьваєет направление силь 
тяжести и сложно для вьтчислений. Если вектор 8 я В любой точке пространства направлен в 


одну и ту же точку Ро ) і -1,3 (вдоль нормали к поверхности зллипсоида), уравнение (1) 
упрощаєтся | Алексидзе, Гелашвили, Картвелишвили, 19721: 


Лк(х,)- 28(х, )/ т з м(х,), ра 





Практические способь перерасчета значений силь тяжести --- постановка и решение гра- 
ничной задачи непосредственно для значений МГПСТ | Алексидзе, 1985|, или решение нели- 
нейной внешней граничной задачи Алексидзе для потенциала силь тяжести: 


яв дді сх, ) хе С 


ДиїЇх. 
ох) 0, хеС 





Линеаризация задачи (3) возможна в предположениий касательно направления МГПСТ на 
поверхности 5, которое вносит значительнье погрешности, а определение направления сильт 
тяжести -- громоздкая задача, как и определениє МГПСТ. Для внешней тяготеющей точки 
получено уравнение | Алексидзе, Гелашвили, Картвелишвили, 1972| 


З дсова, -- 
зАсова, У/8, па обо -0, із-і,3, 
терни Х 


определяющее модуль силь тяжести є по известньїм распределением углов а; в пространст- 


і 


ве С (или косинусьт углов а; по известньм значениям 7? в С с неопределенностью соб а; на 


бесконечности). 

Вопрос разрешимости и единственности задачи (3) для потенциала сильт тяжести принци- 
пиален в гравиметрии: отвечаєт, возможно ли вьтокоточное однозначное определениєе силь 
тяжести во внешнем пространстве по точньм абсолютньм значениям силь тяжести на замк- 
нутой поверхности 5, целиком охватьвающей все массь, и какие нужнь минимальньєе допо- 
лнительнье даннье. Вьчисление силь тяжести во внешнем пространстве при условиях од- 
новременного определения на поверхности 5 ее абсолютньх значений и направления несло- 
жно, но дорого и делаєтся редко. 

Решают задачу (3) путем последовательньх приближений. В плоском случає (граничньтее 
значения (5) и функция и зависят от двух переменньтх) задача определения гармоничной 
функции, модуль градиента которой задан на замкнутом контуре, 





не имеет единственного решения | Алексидзе, 1985|. Если / (2) щ и(х, з їх, ), і-1,2 - про- 
извольная голоморфная в С функция, действительная часть Ке / (2) которой совпадаєт с ре- 
шением граничной задачи (4), ее мнимая часть Пт / (2) также есть решением | Алексидзе, 
1985| задачи (4). Отношение Коши -- Римана К, -К, - - І - - условие перпендикулярности 


ди Іди ду |ду 
градиентов и, у, где К, з 64, 2 - - / - - и К, з ї8а, - - -- / - -- -- угловьгю козффициентьт 
Ох. дх Ох, ОХ, 


касательньгх вдоль градиентов функций и и у (а, а, - - Угль между положительньм направ- 
лением оси и й И у). 


і; ди М 
В конечной уи-окрестности направления / градиента ягад / - З одного из решений гра- 


ничной задачи (4) не существует другого решения задачи (4), т. е. нет двух решений задачи 
(4), направления градиентов которьгх составляли бьт угол меньше цу (угол между градиентами 
двух решений (4) -- 7/2, решения, различающиеся постоянной - - тождественньт). Проверка 
зтого предположения осуществлена путем численного решения методом | Алексидзе, 1978| 
граничной задачи (4) с помощью итераций 
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од ОНИ? 
где козффициентьі а, (и )- . У з направляющих косинусов полученьт нор- 
іх бю 
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мировкой граничньх условий задачи (4), и" -- К-то приближение ее решения. 
Погрешности трансформаций. Погрешность є заменьт уравнения (2) на уравнение Лап- 
ласа Д» (х;) - 0 является разностью решений на сфере С радиуса а граничньх задач для 


уравнений сильт тяжести де(х,)- 22(х фура -0, хє С и Лапласа Лех ) -0 хеС,содина- 





ковьмм граничньм условием 51, - у(5 У Вектор 5 « Направлен в одну и ту же точку О (2, п, 
0). Погрешность удовлетворяет граничной задаче Де - 22 (х,) / г", х є С с однородньми 
граничньми данньми 4, 5 7 0, причем и « є х у. Если начало координат совпадаєт с центром 
сферьт С, то в точке х, «0, і- 18 имеем 


2 2; 
а?у, (64 є ча" у,/б6, (6) 
и при условий, что область С касаєтся поверхности Земли, нормальное гравитационное поле 
будет равно 
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у, з пахов(/т? «тах в(х)/ шіпг| збо 





хєС хеС (6360) 
| 1 2-(978000 --300а) 
РУ 2 ; тах У 
Йаро хієб 28 о ут бі хієЄЄ З (сх У мес 1 й (6380 як а) | 


Отсюда, согласно (6), при а - 100 км получаєм 76 « є х 81 мГал. В перерасчете силь тяжести 
в свободном пространстве по формуле Пуассона погрешность может достичь 100 мГал на бо- 
льших площадях земной поверхности. 


2 2 зі 
В плоском случає у, а?-Уїх 4«єху, а?-Уїх, 4 при х, -0, і - 1,3 дает51 «є х 
ізі ізі 


54 мГал. 
Согласно (6), є 2-2 - 98 2 0, позтому действительное значение сильт тяжести всегда мень- 


ше, чем полученное в предположениий ее гармоничности (на зто указьгваєт и Д? 2 0 - - усло- 
вие вогнутости поверхности). 


Если на сфере С величина 2/ г" постоянная, то погрешность є; являєтся решением грани- 


М 1 
чной задачи Дє, - 28 ечба 


По 


-ургеС,где п) - - среднее значение г в сфере С. Вьтражениє 
я 


1 3 
Є й -Уїх; Ок мажорируєет є;, позтому є; «2,9 мГал, априа х 100 км є, «1 
ізі 


мГал. Таким образом, членьт при 7 в уравнении (2) почти постояннь. Сравнение соответству- 
ющих погрешностей указьгваєт, что уравнение силь тяжести лучше аппроксимировать урав- 
нением (2), чем уравнением Лапласа. 

Уравнение силь тяжести (2) при г - г, бьшстрее сходится, а при малом а/, погрешность 


3 
зтой замень! є, с 22 хз | а" - У, х, | / 5 г) имеет зкстремумь х зх, 20, хз - 2а/ /3. Награ- 


іі 
- Р б з з ни 
нице 5 аномальной области С є» - 0, т.е. тах | Є» | -42а і 15 3 г , Откуда при а - 100 


км є, «0,5 мГал. В зтих оценках погрешности замень уравнения (2) уравнениєм Лапласа 
предположено, что на замкнутой границе 5 известнь точньєе значения силь тяжести; при- 
ближеннье значения 5 х Несколько увеличат величину погрешностей. 

В практике гравиметрии при перерасчете во внешнем пространстве решают не внутрен- 


нюю граничную задачу, а внешнюю (ее сводят к решению внутренних граничньтх задач 
|Алексидзе, 19651). Погрешность перерасчета во внешнем пространстве С (заменьт внешней 


граничной задачи Да(х )- 29 (х у г' з 0 навнешнюю задачу Де, (х, ) - 0 при одинаковьгх 


ста 


1 
граничнькх условиях 8, ш м(5 ) 5), -0)равна є - 5 - 1) и достигаєт максимума в об- 


ласти С при условий РИ шо (а)/ 4, где 2 (а) - - значение силь тяжести на поверхности 


сфероида. Предположение гармоничности силь тяжести в перерасчете на вьтсоту 6000 км 
дает погрешность » 245 Гал. 


0 М М(т- С со80 
Радиальная компонента У, - -- ( 3 ) - анома- 


дл нав на 0) (Р -дсовоч іт) 


лий У(х,) а МИ( хі якхо (ха - 0-7) от однородной сферьт массой М с центром в точке Р (0, 








ПРИ 
0, С « 6375) имеет погрешность перерасчета є, - - - У,, где с -- произвольная постоянная. В 
г 


точках гг аи 0 - 0 погрешность равна Є, - МІ-аКа- 2) г -- У, 1, а ее максимум при любом 
фиксированном 7" существуєт в 0 - 0, возрастая от -оо (7 - 0) до -юо ( хз - 0), спадаєт к отрица- 
тельному максимуму в точке г - г (0) и равен нулю в точке г - а. Если масса М - 100 (а - б)?, 


аномалия в ближней точке поверхности (0, 0, а) равна 100 мГал; приа- 6371,0 - 6361 и га 
6381, є, я 8 мГал. 


Итак, трансформацию наблюденного или нормального поля силь тяжести во внешнем 
пространстве можно осуществлять путем решения методом |Алексидзе, 1978| граничной за- 
дачи для уравнения силь тяжести (2), конкретньг6е операторь которого для разньїх референц- 
поверхностей |Алексидзе, 1985 | в случае земного зллипсоида вращения (с достаточно мальтнм 
сжатием) не генерируют существенно различньх результатов близ Земли. 

Необходимье в численном моделирований явнье аналитические вьтражения фундамен- 
тальньжх решений уравнения силь тяжести (2) полученьт в общем виде |Чернькй, 1982,|. Для 
вьтчислений важно представлять соответствующие интегральт в злементарньх функциях. 

Перерасчет «аномалиий» силь тяжести 02 - є (5) - у (5) (разности нормального и наблю- 
денного поля сильт тяжести) имеет свой нюансь |Дубовенко, 2010). Впрочем, с точностью 0 
« є 2 0,04 мГал аномалиий силь тяжести на равнине является значениями вертикальной сос- 
тавляющей по направлению модуля градиента. 

Перерасчет в нижнее полупространство известной плотности с (Р) в уравненим (2) несло- 
жен. В практике региональной гравиметрии при перерасчете в сторону масс нередко на пове- 
рхности 5 известнь вертикальньгке составнье ОМ /0 хз аномалий, а под поверхностью анома- 
льнье массь расположень на неизвестном расстояний от нее и следует гармонически продо- 
лжить 0У/ / д ха в сторону масс к ближайшим особьгм точкам поля. Задача в подобной поста- 
новке существенно некорректна (имеет бесконечное множество решений | Алексидзе, 19851). 


Результатьт вьічислений. Численно решена |Алексидзе, 1978| нелинейная граничная за- 
дача (3) для пространственного случая с помощью итерационного процесса 








3 ди 
п 2 
Ди, -0, уча и) 509 (5), (7 
тет Ох, 
715 
гдейи, ий, | -- п- и (п - 1)-е приближение, а, (и, 4) - - і-й направляющий косинус градие- 


нта (п - 1)-го приближения и, |. С целью проверки сходимости итерационного процесса (7) 


решень такие внешниєе граничнье задачи: 
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где п а 10, хз - 5 км, 5 - - Горизонтальная плоскость. Точное решение граничньтх задач (8) и 





110) 


(10) таково: 


Пп п п 
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ху ХХ а -р) хх зп) ху 

а граничная задача (9) точного решения не имеет. Рис. І иллюстрирует результатьт решения 

за 5 итераций задачи (3) в 30 точках поверхности 5 на равномерной сетке с шагом Ах; - 10 





иг- 





км в области -25 З х, 225 км, - 25 5 х, 2 15 км (нумерация узлов -- рост х; при нарастаний 
х» ). Первое приближение направления грубое - - а (и) з а, (и,) - 0, дз (и) - 1. В другом 
варианте решения взято а, (и,) - 0,8, а, (и|) - 0, а, (и,) - 0,6, аградиент и, и правая 
часть граничного условия (8) совпали с точностью до є - 10 9заб итераций. Модули точного 
и приближенного градиентов совпадают с точностью 10 7. 

Граничнье задачи (8)- -(10) решали для разньх первьїх приближений и положения вспо- 
могательньтх точек | Алексидзе, 1978), определяющих фундаментальньтг6"е решения. В случає 
рис. 1 вспомогательньте точки взятьт на плоскости хз - 1. Число итераций в вариантах варьи- 


ровало незначительно (не более 10). Граничньге значения точньх и приближенньх градиен- 
тов совпадают с точностью є - 107, они сами совпадают с меньшей точностью из-за плохой 
обусловленности задач с производньми в граничньх условиях (для них, в отличие от задачи 
Дирихле, несправедлив принцип максимума). Зту особенность, не связанную со сходимос- 
тью итераций, следует учитьвать в практическом моделирований. 

В определениий погрешности решения задачи (8) задействовань 36 фундаментальньх ре- 
шений, равномерно расположенньх на плоскости хз - 1, граница 5 - - плоскость ху -5. 
Левая и правая части граничньтх условий (8) после четвертой итерации совпали с точностью 
є з 1079, относительная погрешность решения и - - 10 90, абсолютная --- 10 |; в производньїх 
погрешности достигают 50 Ур. Увеличение числа функций до 49 снизило погрешность 
решения и, производньхх -- до 3-5 906, ди / 0 Хз -- до 20 Ус. Наращивание их числа до 64 


увеличило погрешность решения и производньх до 70 У0, хотя правая и левая части 
граничного условия (3) после восьми итераций совпали с точностью є з 10, Зтот зффект 
обьясняют | Алексидзе, 1985| неминимальностью систем фундаментальньх решений и 
малостью детерминанта матриць козффициентов для фундаментальньх решений. 

Точка х; з Х, - Хз - 0, создающая поле в задаче (8), - - одна из вспомогательньтх |Алек- 


сидзе, 1978), определяющих фундаментальнь6е решения. Зтим обьясняєм немалую точность 
определения производньтх ди / 0 х, по приближенному решению. Если она не входит в состав 
вспомогательньтх точек, точность определения производньх хуже, особенно для ди / 0 хз. 

Если направление градиента на всей границе 5 образуєт с внутренней нормалью острьй 
угол, внешняя граничная задача (3) решается с удовлетворительной точностью и 
погрешность обусловлена не плохой сходимостью итераций (7) (число их мало), а плохой 
обусловленностью граничной задачи с косой производной. 

В общем случає, граничная задача с косой производной вдоль направления І/ 


ди 
--| зм(5) (11) 
5 
при условий, что направление внешней нормали М неортогонально направлению | (І ни в од- 


Ди «0, 


ной точке границь С не совпадаєт с касательной) 
іпб (М, 1)» 0 (12) 
є 5 


х 
не имеет решения, отличного от постоянного |Бицадзе, 1966|. В граничньхх задачах (8)- (10) 
зто условие вьшолняєтся на любой конечной части плоскости 5, так как поле создано точеч- 
ньюм источником, а на всей бесконечной плоскости - - нет, ибо в бесконечно удаленной точке 


плоскости 5 (М, І) | 2. 7 0. Для внутренней граничной задачи (11) и точечного источника 
о о 


Р вне области С условие (12) для достаточно гладкой границьт невьштолнимо, так как всегда 
будет существовать точка на 5, в которой вектор с началом в Р ортогонален к нормали. Чис- 
леннье зкспериментьт на сфере С обнаружили плохую обусловленность граничной задачи 
(11), хотя граничньве условия удовлетворяются хорошо, вьтчисленньт"е и точньєе значения гра- 
диента и и его производньх существенно разньг. 

Решение внутренней граничной задачи (8) для единичного куба С (внешняя нормаль 
имеет шесть разньгх направлений, начало координат (0, 0, 0) расположено на расстояний 0,2 
от нижней кромки куба С, 24 вспомогательнье точки взять на кубе со стороной 1,4) 
согласно (7) при а, (и) з а» (и) -0, аз (и, ) - Ї таково: после пятой итерации граничнье 


значения удовлетвореньт (є « 1 90), но приближение и и его производньве не имеют общего с 
точньм решением (/ « в/ хі ху ху иего производньми; увеличение числа итераций до 


24 не увеличило точность приближений. В табл. 1 в четьтрех точках границьт (точки 
коллокации) даньт точньєе и приближенньте значения градиента и производньтх (в других 
точках они из-за симметрии совпадают с одним из приведенньтх в табл. 1). Добавив 25-е 
фундаментальное решение (оно совпало с точньм решениєем задачи (8)), после четвертой 
итерации граничньєе условия, решение и производньєе мтали удовлетворять с точностью є - 
107. Чтобьг разделить влияние неустойчивости граничной задачи (11) на суммарную 


погрешность решения граничной задачи (3), для функции Ї/ / хі ху ху решена 


граничная задача (11). Последние столбць табл. 1 дают решение и граничной задачи (11) в 
тех же точках: погрешность обусловлена неустойчивостью задачи (11), а не медленной 
сходимостью итераций (7). 

Сходимость итераций (5) в плоском случає имеет особенности: точность решения 





внешней граничной задачи (4) с граничньгтм условием у(5) з у Ух хо | є1 иточньм 
5 


решением и - іп 4 хо "б х; для единичного круга с центром в начале координат при на- 
чальном направлений а? - со8 (г - ф). ат - 8іп (1 -- Ф) (1-- центральньй угол, ф-- 
возмущение начального направления) зависит от Ф и М. В табл. 2 приведено количество 
итераций М до условия их остановки | Є | з. | | ягад и | - | « 107? , максимальнье 


погрешности решений в точках коллокации после каждой итерации для производньхх - - 
є « 1077, а между точками коллокации в конечном решений є, я 10, где 


1,2 
ди ди? ди 
єї | | -|етах| оо 


Ох, дх. дх. 


Решение внутренней граничной задачи и при благоприятньтх условиях определяется с 
точностью до постоянной (позтому не сравниваєм точное и приближенное решения) и знака 
(при Ф - л первая итерация удовлетворяєт граничной задаче). 


Решена задача (4) с радиусом круга г - 2 и направлением градиента а - є - /2 /2 


(совпал с прямой х; - х, ). Зто направление образуєт с нормалью угол 907? в точках круга 


Ма) ооо ри М, | ото, то |, т. е. не вьшшолняется условие (12). Точное решение 


задачи и - 0,693, а граничное условиє м (5) - 0,5. В приближенном решений задачи (4) (рис. 
2 - - результат девяти итераций, где х,, Х, - - кООрдинатьт граничньтх точек коллокации) гра- 
диенть приближенного и точного решения близки, а сами решения и их производнье сущес- 
твенно разньге, так как приближениє а - а9 - 4/ 2/2 не обеспечиваєт условия (12). 


Для сходимости итераций (5) необходимо не только, чтобьі точное решение удовлетво- 
ряло условию (12), но и вьібор начальньїх приближений, для которьх в итерациях (5) не во- 
зникали бь козффициенть а! , не обеспечивающие условия (12). 


Решение по итерациям (5) внутренней граничной задачи для единичного круга с центром 





в начале координат, граничньтм условием м (5) - 1 / У (х - 1) Рв (х» - 1) ? при направле- 


ний градиента аб - 0,8315, а? - 0,5556 (как любоге постоянног направлениєе, зто не обеспечи- 
ваєт вьшолнения условия (12)) иллюстрируєет рис. 3 (М - 8 итераций). Точное решение 





и - т (а - 1) еВ (х» - 1) би комплексно-сопряженная (а з и - іу- - голоморфная функ- 
ция) с ним функция у по соотношению Коши - - Римана также удовлетворяєт зтой задаче. 
При о? итерации (5) стремятся не к точному решению и, а к комплексно-сопряженной функ- 
ЦИИ У. 

Вьтчисление аномалий. На практике отьшскивают не потенциал у/ - и 4 у силь тяжести, а 
потенциал и аномалиий силь тяжести (потенциал нормальной силь тяжести у известен). В ли- 
нейньх приближениях граничньгх задач грави- и магнитометрии переход к аномалий не 
изменяєт не только уравнения, но и вид оператора граничньх условий. Необходимо лишь в 
качестве граничной функции м (5) взять аномалию силь тяжести. 

В нелинейньх граничньх задачах переход к аномалиий изменяєет левую часть граничнькх 
условий (3), (4), так как операция вьтчисления модуля градиента не дистрибутивна: 
| рай у | же 





ягад и | . | ягай у | . Если граничная функция м (5) - - модуль градиента 


потенциала у/ сильт тяжести, для потенциала аномалиий и получаєм нелинейную граничную 
задачу: 
(ди) (оду д 
и у ди 2 2 
Ди «0, У |--| «25 |-- --| зу (5)- |та у| , и, з у,, (13) 
єп |дх, еаногій  2 Я я 
граничньєе значения которой при нулевом потенциале нормальной силь тяжести совпадают с 
граничньми условиями (3). Решают зту задачу одним из последовательньїх приближений 
(условия на оо опущень(т): 





3 0 
У 2 ду і ди, Ин 











еф Оу 0 зо | вгад | 
Ди, 70, зум (5)- -----ччое, 
| ягад (у чи,)| у (5) 
5 
53 Й ду М ди, | дич | 
прі ЧИ 7 З С 2 ЛВС | ягай у| 
Ди, і 2-0, - у (5)- , (14) 
у (5) у (5) 








3 
п 2 2 
Ди, 2 0, У; 12 - - у" (5) - | ягад у| 
не дх,  дх, дх; 
Определение численной скорости сходимости итераций (14) имеет практический интерес, 
так как на каждой итерации решаєтся граничная задача и слабосходящиеся процессьт сущест- 
венно ухудшают зксплуатацию алгоритмов. Сходимость (14) вьічислена на модели 


Пп 1 





Ше У , у (5) є | ага (/ -- У)| на границе 5 
2 2 2 По 2 2 
ху ВХ; ВХ «ОКХ з (з кА) 
(плоскость хз - 5 на 36 узлах коллокации в точках, х; - -25 4 10К, х, « -25 4 10), ха с5,К,) 


- 1,2, 3, 4, 5) для 36 фундаментальньх решений, полюса которьгх находятся под узлами кол- 
локации на плоскости хз - 1) до вьшолнения условия плах | | ягад (у чи,) | - м (5) | «є. При 
5 


є-107 первьй итерационнь6й процесс завершился за три итерации, второй --- 36, последний 
-4 
- 103 итерации (при є - 10 -- за 4, 78 и 146 соответственно). 
Для всех итераций (14) первой решена граничная задача с производной вдоль хз на грани- 


це. Вспомогательньтєе точки для фундаментальньх решений взятьт на плоскости х. - 1. При- 


ближениєе решения и его производньгх из всех итераций хорошо совпадают между собой (при 
є 2 10 сточностью 10, приєс- 10" с точностью 107), но далеки от точньх. Их средне- 
квадратичнье значения в узлах коллокации (1-й столбец) и в середине между ними (2-Й сто- 
лбец) даєт табл. 3, а соответствующие погрешности для разньх є - то» (к-1,2,..3)миМ-- 
табл. 4. 

Последняя таблица показьваєт, что увеличение точности аппроксимации граничньіх ус- 
ловий в узловьїх точках коллокации без увеличения количества узловьгх точек (увеличения 
числа функций в разложениий) не влияєт на погрешность (начиная с К - 2). Зто говорит о не- 
линейном поведений погрешностей решения граничньтх задач с производньми и побуждаєт 
исследовать погрешности решений каждой конкретной задачи. 

Вьводь. Трансформацию значений поля силь тяжести во внешнем пространстве можно 
осуществить путем решения граничной задачи (3) для уравнения силь тяжести (2), конкрет- 
ньй вид которого зависит от избранной модели Земли. В плоском случає зта операция сво- 
дится к решению граничной задачи (4) по итерациям (5). Точность ее решения существенно 
зависит от мерьт обусловленности задачи, приближений направления искомого градиента, 
параметров численного метода | Алексидзе, 1978| (числа и положения фундаментальньіх ре- 
шений), направления внешней нормали, вдоль которой вьтчисляют производнье в граничньтх 
условиях (вьшолнение условия (12)). Невзирая на зто, ее решение на простьїх моделях при 
определенньх ограничениях успешно. Основньг6е источники неоднозначности решений 
задачи - - конечное вертикальное разрешение решений, неадекватньгй вьтбор модели задачи и 
начальньх приближений решения. 

Благодарю академика В. ИЙ. Старостенко за ценнье замечания касательно рукописи ста- 
тьми. 
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Рис. 1. Решение внешней граничной задачи Алексидзе (3) и фрагмент соответствующей 
таблиць данньтх (х 107). 

Рис. 2. Решениє внешней граничной задачи (4) для кругасг «2 при ад з ад - 42 /2. 


Рис. 3. Решение внутренней граничной задачи (4) для единичного круга. 


Таблица 1. Решение внутренней граничной задачи (8) и задачи с косой производной (11) 
для единичного куба 


Но- и ди ди ди ии ди, | ди, | ди, и ди ди ди 
Че 0 з 0 - 0 з дх д х. дха 0 0 - д хз 











ромюю 





За 39 4 4 45 б 
і 28 0 0 8 ї 

і 36 18 2 8 6 й 
30 18 0 8 2 


Таблица 2. Параметрь решения плоской внешней граничной задачи (4) 
для единичного круга 





Таблица 3. Приближеннье решения нелинейной граничной задачи (20) 
на плоскости 


ди/дхя 
0565 | 098 | 0,029 | одрог | 0029 | 0021 | 0026 | 0067 
0605 | обії | 0032 | 0028 | о0з2 | 0028 | 0038 | 0053 
0617 | 0614 | 0030 | 0032 | оозо | о0з2 | обат | 0,43 





Таблица 4. Погрешности решений нелинейной граничной задачи (14) 
на плоскости при разной точности 


є(диідхі) є(ди/дхі) є (ди/дхз) 2 





Аналитическое продолжение сильт тяжести в региональном масштабе определяет решениє 
нелинейной граничной задачи Алексидзе для уравнения силь тяжести. Дана оценка погреш- 
ностей ее заменьг на соответствующую задачу для уравнения Лапласа. Точность ее решения 
существенно зависит от мерьт обусловленности задачи, приближений направления искомого 
градиента, параметров численного метода, направления внешней нормали, вдоль которой 
вьічисляют производнье в граничньх условиях. 


